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Abstract 

We study the transmission through a double quantum-dot system in the Kondo regime. An exact expression for the 
transmission coefficient in terms of fully interacting many-body Green's functions is obtained. By mapping the system into 
an effective Anderson impurity model, one can determine the transmission using numerical renormalization-group methods. 
The transmission exhibits signatures of the different Kondo regimes of the effective model, including an unusual Kondo 
phase with split peaks in the spectral function, as well as a pseudogapped regime exhibiting a quantum critical transition 
between Kondo and unscreened phases. 
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1. Introduction 

Semiconductor quantum dots have played a major 
role in the investigation of strongly correlated effects 
in nanoscale systems, as highlighted in the pioneering 
experiments on Kondo effect in these devices [1] . More 
recently, double-dot setups have been used to investi- 
gate two-impurity and two-channel Kondo physics [2]. 
These experimental developments have spurred great 
theoretical interest in double dots in the Kondo regime. 
In particular, studies of side-coupled [3] and parallel 
[4] dot configurations have been carried out. 

Motivated by these accomplishments, we study 
transport properties of double quantum dots (DQDs) 
with one interacting dot ("dot 1") coupled to a large 
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dot, effectively noninteracting ("dot 2", with energy 
£2). This seemingly simple setup can be described 
by an Anderson model coupled to a fermionic bath 
with a nonconstant density of states (DoS). Compared 
with the standard case of a metallic host (constant 
DoS), a nonconstant DoS leads to nontrivial new and 
interesting features in the many-body Kondo state. 

As shown previously [5] , a splitting the Kondo res- 
onance appears when the DoS shows a sharp peak at 
the Fermi energy, while the Kondo singlet itself is pre- 
served. The proposed DQD setup also allows for the 
appearance of a pseudogap in the effective DoS, lead- 
ing to a critical transition between Kondo and non- 
Kondo phases. These phenomena substantially mod- 
ify the spectral function of the interacting dot. We 
show here that another quantity, the transmission co- 
efficient, can also be used to explore the critical tran- 
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sition and other features of this system. 

We obtain an exact expression for the transmission 
coefficient T{uj) in terms of the fully interacting many- 
body Green's functions for the DQD system, which are 
then calculated within the numerical renormalization- 
group (NRG) framework. This approach allows us to 
fully describe the characteristics of the transmission 
coefficient in the different regimes mentioned above. 



2. Model 

The model describes the DQD setup depicted in Fig. 
1: quantum dots 1 (singly occupied with charging en- 
ergy U and energy ei) and 2 (energy £2) are coupled 
to metallic leads {£ = L, R) and to each other. 




Fig. 1. Schematic representation of the double-dot system. 

The Hamiltonian is given hy H = Hi + H2 + -ff Leads + 
-ffoot-Dot + -ffoot-Loads, where 

k, l=L,R 

ii"i=i,2 = Eicl^c^^ + 5i,iUni-[n\i , 

-^Dot-Leads = ^ iyUi'-'^laCfka + "^^Wka 4k<T ) , 
k, l=L,R 

//Dot-Dot = Ac|^C2<, + A*4^Ci^ . (1) 

The transmission amplitude for an electron tunnel- 
ing from the left lead to the right lead can be ob- 
tained from the retarded Green's function connecting 
the leads, G'Lkflk'(^) = {{^l^, c^k')): 

GLkflk'(t^) = ^ ffik V'ilkGij (Lj)y;flk' 5Hk' • (2) 

ij = l,2 

with gosL = {^^^ — £ik)^^ ■ In the derivation of the ex- 
pression above, we have used identities from the equa- 
tions of motion for fermionic operators Ca, in the 
frequency domain that include [6]: 



u^{{Ca:cl))^ = {[c,,cl\,) + {{[c^,H\-:cl)) , 

^((Ca :4))„ = {[c„,ct],)-({c, : [ct,i/]_)). (3) 

In the wide-band limit D ^ \u)\, where D is the half- 
bandwidth, we approximate Vifk = Vu (k-independent 
couplings) and X]k('^^ ~ ^fc)"^ ~^ —ii^po (where the 
same DoS po is assumed for each lead) and one ob- 
tains the following expression for the energy-dependent 
transmission from L to R: 

T{w) = Gl^^rw (^) = 27vpoJ2 V*LG,jiu:)VjR . (4) 

k,k' ij = l,2 

In the following, we assume a symmetric configura- 
tion ViR = ViL ~ Vi/ \p2 with real values for the cou- 
plings for simplicity. Defining A; = -KpoV^ and A12 = 
TrpoViVb, T(lo) can be written in a compact form: 

T(i^) = l\xGxY A2G22 + A12 (G12 + G21) . (5) 

All Gij in (5) are fully interacting Green's func- 
tions. Interaction effects are introduced into the dot-2 
Green's function by direct and indirect (via the leads) 
tunneling to dot 1. The identities (3) can be used to 
establish the relations 

G22(t^) = + G2i(c^)Gn((^)Gi2(a;) , 

Gi2(a;) = G<°' (a;) (A - iAi2) Gn (a;) . (6) 

Notice that for real couplings, Gi2(w) = G2i(c<j). These 
relationships in turn lead to an exact expression for 
the transmission involving only the noninteracting dot- 
2 Green's function 6*22'' (<^) = (a; — £2 + jA2)~^ and 
the fully interacting dot-1 Green's function 6*11(0;) = 
((Ci.:cL)): 

T{uj) = AiGii(a;) + 2A12 [G§2(t<j)(A-iAi2)Gii(6^)] 
+ A2 G§2(^) [1 + G§2(^^)(A - jAi2)'Gii(a;)] . (7) 

The dot-1 Green's function Gn is obtained from 
NRG calculations in the following manner: Hamilto- 
nian (1) is mapped [5] onto an Anderson impurity con- 
nected to an effective nonconstant DoS through a hy- 
bridization function 

A((jj) = np2{w) 

with p2{oj) = A2/{7r[(a; — £2)^ + A^]}. In general, 
A(0) 7^ gives the effective broadening of the single- 
particle Hubbard peaks in dot 1 produced by the cou- 
pling to the effective nonconstant DoS. Notice that, for 
A2, A ^ (single-dot case), A(0) Ai. 
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Fig. 2. (color online) Spectral densities (left panels) and 
transmission amplitude (right panel) for the side dot con- 
figuration with £i = —U/2 and A = 0.04_D (dashed lines), 
0.06-D (•) and 0.08-D (thick solid line). Notice the appear- 
ance of secondary peaks in |T(lj)P as A increases. 

We focus on two limiting configurations. In the "side- 
dot" limit (Ai = 0; A2,A 7^ 0), dot 1 is connected to 
the leads only by second-order tunneling processes me- 
diated by dot 2, and A(aj) has a Lorentzian form, with 
a peak of width A2 centered at £2. In the "parallel-dot" 
limit (Ai,A2 7^ 0; a = 0), dot 1 and dot 2 are con- 
nected only by indirect tunneling through the leads. 
Now A(ti;) has an "inverted Lorentzian" shape, vanish- 
ing as loj — e2p at = £2; if dot 2 is in resonance with 
the leads (£2 = 0, the case assumed throughout the re- 
mainder of the paper), this corresponds to a pseudogap 
in the effective DoS with exponent 2. 

The spectral function A\\[uj) = — 7r~^Im Gii(ti;) 
is obtained from the NRG spectra using the method 
described in Ref. [7]. In short, at each NRG itera- 
tion A'', one collects partial information for the spec- 
tral function by approximating Ai\{uj) — y4f^(aj) for 
UJN < < lOtJjv where ujn ~ DA~^^^ (A > 1) is 
the characteristic energy scale at iteration A^. A con- 
tinuous function can be obtained by replacing the 5- 
functions at energies En entering the Lehmann expres- 
sion by logarithmically broadened lines S^{lo — En) oc 
exp [-(Inuj -InEnf/b^]. 

Next, Re Gn (lo) is obtained from An (lo) via the ap- 
propriate Kramers-Kronig (KK) transformation. Be- 
cause of the subtleties involved in the use of such trans- 
formations, we compared (i) the results of transforming 
Aii{ll}) obtained numerically as described above with 
(ii) those from performing the KK at each NRG step. 



calculating Fn{uj) = n^'^'P J dio' Aii{lu')/ {lu' - u) 
and then approximating ReGii(tj) — — 7r^^i<jv(aj) 
for uiN < ^ ^ lOtJjv. We obtained qualitatively 
similar results with both methods, with better agree- 
ment for larger b entering 5''{u} — En) {b — 0.3-0.6 
for the results presented). For large A'', an addi- 
tional, broadening-independent estimate of Gii(O) 
can be obtained from the expression [8] Fn{0) = 
^^(0)-' E„ {\{n\ci.\0)\^ - KO|ci.|n)n /E^, where 
En , \n)N and Zn{0) are, respectively, the eigenener- 
gies, eigenstates and zero-temperature grand-canonical 
partition function at NRG iteration A', with n = 
representing the ground state. 



3. Results 

In the following, we take U ^ 0.5D and A2 = 0.02D. 
In the side-dot configuration considered, the effective 
dot-1 hybridization A{u)) has a resonance of width A2 
at the Fermi energy. Note that A(0) = A^/A2, so that 
the Kondo temperature Tk oc e^^'''^^"-' will increase as 
A increases. 

As previously shown [5] , the presence of a resonance 
in the effective hybridization leads to a splitting in the 
spectral function for A > A*, where A* is defined im- 
plicitly by Tk{\* , A2) = A2. This can be seen in the 
left panels of Fig. 2, for which A* ~ 0.05D. 

The transmission amplitude is shown in the main 
panel of Fig. 2. At = 0, the Friedel sum rule requires 
that ImGii(O) — — A2/A^. Substituting this result in 
Eq. (7) leads to |T(0)P = independently of \, A2. 
The vanishing of transmission across the double-dot 
system is essentially an interference effect: the path go- 
ing directly from the left lead to the right lead through 
dot 2 interferes destructively with the alternate path 
in which the electron tunnels in and out of dot 1. 

Signatures of the splitting in the spectral density 
appear in the transmission amplitude at finite ui. For 
increasing values of A, the onset of the splitting in the 
spectral function (left panels in Fig. 2) is accompanied 
by the appearance of satellite peaks in \T{lu)\^ near 

LU = 0. 

In the case of the dots in a parallel configuration, 
the situation is different, as shown in Fig. 3. Notice 
that |r(0)p — 1 always, which is a direct consequence 
of the presence of the pseudogap in the hybridization 
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Fig. 3. (color online) Transmission amplitude for the par- 
allel configuration with Ai = 0.05D, and ei = —U/2 
(dashed line) -0.088C/ (•) and -0.086C/. The critical point 
is reached at = —0.086524(7. Inset: A dip in T{uj) crosses 
a; = due to the pseudogap-induced transition. 

function A{uj). In the resulting effective pseudogapped 
Anderson model, the spectral density vanishes as 
at the Fermi energy. Thus, Im Gii(O) = 0, effectively 
closing the transmission through dot 1 at the Fermi 
energy. In this regime, the transmission is dominated 
by the resonant tunneling through dot 2. 

Important differences appear in particle-hole (p-h) 
symmetric and asymmetric regimes. In the p-h sym- 
metric case (ei = —U/2, dashed lines in Fig. 3) the 
dot-1 spectral density varies as lo'^ within a relatively 
large range < A2 around the Fermi energy [5]. In 
this case, |r(6t;)p is essentially a Lorentzian of width 
A2. 

Away from p-h symmetry, a quantum critical point 
separating Kondo and non-Kondo phases can be 
reached [5]. Passage through the quantum critical 
point is reflected in position of a peak in the dot-1 
spectral function at a frequency u* that crosses from 
a;* > in the Kondo phase to tj* < in the unscreened 
phase [9]. 

We find that the transition also has a signature in 
the transmission. The inset to Fig. 3 shows that the 
low-energy peak in An translates into a dip in jT(tj)p 
at (j; « u)* — a dip that passes through the Fermi en- 
ergy at the quantum critical point. The amplitude and 
width of this dip depend on structural parameters that 
should be tunable in experiments to enhance the fea- 
ture. This opens interesting possibilities for controlled 
experimental study of a quantum phase transition 



through conductance measurements. 



4. Conclusions 

We have analyzed the electronic transmission in two 
different regimes of a double quantum-dot system. Uti- 
lizing an Anderson Hamiltonian and its exact solu- 
tion using numerical renormalization group methods, 
one can determine the energy-dependent transmission 
function for the structure. As different geometries are 
explored, one can access an unusual Kondo regime with 
split peaks in the spectral function, as well as a Kondo 
system in a pseudogapped environment, allowing ex- 
ploration of an interesting quantum critical transition. 
The transmission function exhibits clear signatures of 
these different Kondo regimes, opening the possibility 
of extensive experimental studies of their properties 
and response to different perturbations. 
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